A study on bifurcation diagrams in relation to synchronisation in chaotic systems 
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We numerically study some of the 3-D dynamical systems which exhibit complete synchronisation 
as well as generalised synchronisation (GS) to show that these systems can be conveniently parti- 
tioned into equivalent classes facilitating the study of bifurcation diagrams (BDs) within each class. 
We demonstrate how BDs may be helpful in assessing the robustness of GS and in predicting the 
nature of the driven system by knowing the BD of driving system and vice versa. We extend the 
study to include the possible GS between elements of two different equivalent classes by taking the 
example of the Rossler-driven-Lorenz-system. 
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I. INTRODUCTION 



Ultra-sensitivity on the initial conditions makes it hard 
to predict the future behaviour of a chaotic system which 
is termed to have long term unpredictable behaviour. 
The question if two chaotic system running side by side 
may be forced to follow the same (or functionally re- 
lated) paths on the respective attractors did raise inter- 
est in the scientific community long back with the works 
of Yamada and FujisakafEi; and Afraimovich et. oi.0 
who dealt with coupled chaotic systems. The subject 
of synchronisation was thus brought forth. The theory 
of synchronisation was given a polished look by Pecora 
and Carroll [J] who along with the seminal paper by Otts 
et. al. Q have simulated wide range of research activities 
that has spanned a period of almost two decades. With 
possible applications in the fields of electronics, secured 
communication system, lasers, chemical reactions, biol- 
ogy, robotics etc., synchronisation is, no doubt, a topic 
of recent interest. For details see review by Pecora et. 
aL0] and the references therein. Various possible types 
of synchronisation are: complete, lag, phase, burst and 
generalised (GS); the last one was introduced for unidi- 
rectionally coupled systems in a paper Q by Rulkov et. 
al. Kocarev and Parlitz[8| went on to develop a general 
theory for GS and showed that it implied predictability. 
Also, mutual exclusiveness of GS and conditional equiv- 
alence was showcased. The paper also touched upon the 
influence of parameter mismatch on GS and hence, ro- 
bustness of GS was discussed. 

Our interest here is to introduce a systematic 
study of the relationship between bifurcation dia- 
gram (BD) and GS in unidirectional coupled sys- 
tems (master (drive/driving)-slave (response/driven) sys- 
tems) in the literature of synchronisation and to discuss 
its important applications, e.g., assessing the robustness 
of GS etc. Numerical experimental results will also be 
presented in due course to validate the claims made in 



the paper. 



II. EXTENDED SYNCRONISATION 



Let the following two equations be the driving system 
and the driven system respectively: 



x = f(x;p x ) 

y = g(y,h(x);p y ) 



(1) 
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where, x £ R™, y £ R m , p x and p v denote the cor- 
responding sets of parameters, and h(x) is the func- 
tion responsible for driving the driven system. Now, 
for a given value of p x and p y GS is said to exist for 
(JTJ) if 3 a transformation H Pi . P]j : R™ — > M. m , a man- 
ifold M — {(x, y) 3 y = H PtiP (x)} and a subset 
B = B x x By C W 1 x W n with M C B 3 all the tra- 
jectories of ([1]) with initial conditions (xo,yo) in basin B 
approach M as t — ► oo. Obviously, complete or identical 
synchronisation (IS) means H PxtP is an identity trans- 
formation. The necessary and sufficient condition^ for 
GS to occur is that V(x ,yo) S B, y = g(y, h(x);p y ) is 
asymptotically stable. This means one can predict the 
behaviour of y, based on the knowledge of x and ~H Pj .. p ■ 
Also, one may define the vector field y to be condition- 
ally equivalent to x if 3 a C k diffeomorphism which takes 
orbits of f to orbits of g, preserving the senses. Again, 
suppose that IS occurs for p x — p y = po (say), and if 
3 a neighbourhood P of pq 3 \fp y € P driven system 
is asymptotically stable, one would say that GS persists 
and hence GS is robust. 

Now we, for the sake of convenience without losing 
generality, fix all the parameters of driving system and 
driven system except one for each and let p x and p y de- 
note these two. We keep p y fixed and vary p x continu- 
ously and suppose that GS persists, unlike IS, for p x 6 P. 
This independence of GS from p x implies some sort of ro- 
bustness which in general need not be same as the one 
define earlier. To be more concrete, we illustrate our idea 
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with the following Lorenz system as the driving system: 
X l = <t(x2 - X\) 

x-2 = rx\ — X2 — X1X3 (2) 
x 3 = xix 2 - bx 3 

which drives a similar driven system with y\ replaced by 
x\ and r replaced by r — Ar. If one defines e = x— y, then 
the evolution of the trajectories in the error space defined 
by e would depend on Ar. This means if for p x G P, such 
a driving system keeps the driven system in generalised 
synchrony that becomes IS when p x — p y = po G P, then 
effectively we have a family of systems (same coupled dif- 
ferential equations but different parameter value of a cer- 
tain parameter) such that any two can be in GS. This can 
be checked by plotting, for the entire P, the maximum 
conditional Lyapunov exponent which of course should 
be negative. This can be taken as an alternatively as 
the definition of P. One such plot (fig. (1)) has been 
illustrated for the Rossler system. There is no dearth 
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FIG. 1: ES in unidirectionally coupled Rossler system: max- 
imum conditional Lyapunov exponent has been plotted with 
the parameter b. (One may refer fiol ] for the standard form 
of the system.) 

of such families as our numerical experiments with sys- 
tems like Lorenz, Rossler, Chen, Rucklidge, Nose-Hoover, 
Rabinovich-Fabrikant equation, Duffing's two well oscil- 
lator etc. would suggest. For such a system, we shall say 
that the system possess extended synchronisation (ES) 
with respect to the parameter p (say) . 

Let us call such a particular family F p (where the su- 
perscript stands for the neighbourhood of the parameter) 
and its i th element ft; one may say fi and fj are extend- 
edly synchronised to each other. One may note that the 
relation (say, E)-"is extendedly synchronised to"- is an 
equivalence relation, for, it satisfies 

1. Reflexivity: fiEfi V/, G F p . 

2. Symmetry: /,/•./, > J, I: J, V/,;,/, G F p . 

3. Transitivity: fiEfj and fjEfk 



What this amounts to is also mathematically 
interesting^. It says that if we form a set <f> of 
all such families F Pi then each F Pi is an equivalence 
class and forms quotient space. We shall call such 
a family extendedly- synchronised- family (ESF). Any 
system fi in a particular ESF, F i 1 C $ would then be 
the representative of the ESF. 



III. GS, ES AND BD 

Now suppose we pick a particular ESF, F (say) and 
choose arbitrarily an element of it f driven (say). The 
question: can we induce bifurcation in the system f driven 
without changing its parameter which is kept fixed at po? 
The answer is, yes. We just need to couple f driven with 
any other element of F p , f driving (say) and continuously 
vary the parameter, p of f driving 3 P £ Pand then plot 
a variable of f driving vs. p and the corresponding conju- 
gate variable of f driven vs. p to find that both the BDs 
are almost exactly identical. See fig- (2) and fig. (3) and 
their captions for examples validating this point. The 
more important and interesting point to note is that one 
has, thus, induced bifurcation in f driven sans tampering 
with its parameter pq. This striking alikeness of BDs 
may be used as a very crude pictorial definition of ex- 
tended synchronisation. Also, it may serve to illustrate 
the robustness of GS, a property which requires to find 
a Lyapunov function, not an easy thing to do in general. 

We now ask a deeper question: mathematically, when 
should two BDs be alike for the cases we are deal- 
ing with? To answer this question let us first note 
that identical BDs would mean the existence of a fam- 
ily D p of C k diffeomorphism, each element of which is 
dp. Pa 3 D p = {d PtPo 3 p G P} because the numerical ex- 
periments show that period changing in the driving sys- 
tem is causing same period changing in the driven system 
which, of course, hints that within an ESF (at least for 
those which have been experimented on numerically) not 
only GS but also conditional equivalence prevails (which 
can be checked by plotting the conjugate variables for 
an ESF for any parameter mismatch such that the pa- 
rameters remain within P and observing that there is 
no 'kink' in the plot confirming at least C 1 diffeomor- 
phism; an example will follow shortly). Then, similar 
BDs would mean that for a particular Poincare section 
in the driving system 3 at least one such Poincare sec- 
tion in the driven system 3 3 a, family, M p of homeo- 
morphisms m PtPa 3 M p — {m PtPo 3p£ P} between the 
two Poincare sections. This introduces predictability in 
a broader sense: knowing BD of driving system one may 
predict the BD of driven system and the maps being bi- 
jective facilitates the existence of corresponding inverse 
maps allowing one to predict the BD of driving system 
from the knowledge of the BD of driven system w.r.t. 
parameter of the driving system. 

As an aside, we would like to highlight a very inter- 
esting coincidence which may be pregnant with fruitful 
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FIG. 2: Alikeness of BDs: figs, (a) and (b) show respectively 
the BDs of the driving and the driven Lorenz systems in com- 
plete replacement (CR) of xi coordinate; r is varied for the 
driving and for driven it is kept fixed at r = 28. a and b are 
kept fixed at 10 and 8/3 respectively for the entire ESF. Sim- 
ilarly, in figs, (c) and (d) BDs are for Rossler systems with 
unidirectional coupling in X\ coordinate with coupling con- 
stant 1.6; b of driving is varied, a and c are kept fixed at 0.2 
and 5.7 respectively. Likewise, BDs of drivings and drivens 
of Chen system (figs, (e) and (f); CR in X2) and Rabinovich- 
Fabrikant equation (figs, (g) and (h); CR in X2) have been 
plotted. For all the figures, parameter is along the abscissa 
and variable is along the ordinate. For the standard forms of 
the equations one may refer 



research possibilities. In all the examples of ESF, on 
which numerical experiments have been done by us, the 
corresponding representative system may be written as 
a jerk eauation\l(\. (One may recall in the passing that 
in general a set of three coupled first order differential 
equation can't always be reduced to a jerk equation.) 
e.g., Lorenz: x\ + (1 + a + b — ii/xi)xi + [6(1 + a + 
x\) — (1 + a)xi/xi\xi — ba(r — 1 — x^)xi = 0, Rossler 
(with a = b — e): x 2 + (c — e + ec — x 2 )x 2 + [1 — ec — 
(1 + e 2 )x 2 + ex 2 ]x 2 + {sx 2 + c)x 2 + e — 0, Rucklidge: 




3.2 3.4 3.6 3.8 4 4.2 4.4 4.6 4.8 5 




FIG. 3: Alikeness of BDs (continued): Two more examples 
of BDs of drivings and drivens of Rucklidge system (figs, (a) 
and (b); CR in X2) and Duffing two well oscillator (figs, (c) 
and (d) CR in X2). The striking alikeness is obvious. In the 
case of Rucklidge system the only parameter of the equation 
is varied. While, Duffing two well oscillator's BDs are w.r.t 
the variation of the angular velocity. For the standard forms 
of the equations one may refer • 



x-2 + (1 + K — x 2 /x 2 )x 2 + k(1 - x 2 /x 2 ) + x\ = 0, Nose- 
Hoover: x'i — (xi + x\)xi/xi — (a — 1 — x\)x\ = etc. 
This means only one variable alongwith its derivatives is 
enough to carry all the information regarding the change 
in parameter being done to construct BD and unidirec- 
tional coupling of the driving system's "single" variable 
with the driven system is enough to induce the entire 
effect of the evolution of the driving system with the pa- 
rameter into the driven system's BD w.r.t. the parameter 
of the driving system so that the two BDs are alike. 



IV. BD IN ROSSLER-DRIVEN-LORENZ 
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FIG. 4: Mimicry of Rossler system's periodic multiplications 
by Lorenz system: figs, (a) , (c) and (e) are of Rossler system 
and corresponding figs, (b), (d) and (f) are of Lorenz system. 

There remains one interesting possibility. It can very 

p p 

well happen that given two ESFs, F i 1 and F- 3 the rep- 
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FIG. 5: BDs of Rossler-driven-Lorenz: figs, (a) and (b) are 
for k — 40 and are showing how BD for Rossler (fig. (a)) is 
causing BD for Lorenz (fig. (b)) to be alike, while for k = 2 
one may see that BD for Lorenz (fig. (d)) is not following the 
BD for Rossler (fig. (c)). and j/2 have been plotted w.r.t. 
64. 
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This arrangement is known to show GS for a = b = 
0.2, c = 9.0,cr = 10, r = 60, & = 8/3 and most impor- 
tantly k = 40. In the passing, it may be noted that 
the driving system may again be written as a jerk equa- 
tion as mentioned earlier. To capture this GS, one may 
painstakingly go on to calculate the continuity and differ- 
entiability statistics [ll| developed by Pecora et. al. Of 
course, the process is rewarding with computers making 
the job easy. With the parameter values fixed at the val- 
ues mentioned above, one may note (as in fig. (4)) how the 
period doubling and the period tripling in Rossler system 
brought about with the variation in c is mimicked by the 
Lorenz system whose parameters are kept fixed. Also, 
the time series of the driving system in non-chaotic state 
stops the chaos in the driven just as the time series of 
the driving system in chaotic state can initiate the chaos 
in the driven even if it initially was in non-chaotic state. 
This mimicry is gradually lost when k is reduced as is 
evident from the fig. (5). Hence, even in this case BDs 
have been able to showcase the existence GS when it re- 
ally exists and this is in consistence with the statistical 
technique used by Pecora et. al. As with the GS between 
the elements of a single ESF, the existence of conditional 
equivalence may also be claimed here because the maps 
for various values of b are smooth. At typical values of 
b (fig. (6)), one may note that there are no kinks in the 
plot ensuring at least C 1 map. 




FIG. 6: Smoothness of the maps in Rossler-driven-Lorenz for 
k = 40: figs, (a) and (c) are respectively the plots for y2 — X2 
for the value b — 0.2 and 0.75 where there is chaos and period 
four respectively in the BD of Rossler (fig. (5)). The figs, 
(b) and (d) are the blown up image of the boxed regions in 
the figs, (a) and (c) respectively, showing that there are no 
'kinks'. One should observe that the plot is not a line-plot 
but a dot-plot to appreciate the fact. 

P P 
resentative elements £ i 7 ^ 1 and fj £ Fj 3 respectively 

may depict GS when coupled. For example, we can have 
a Rossler system driving a Lorenz system through diffu- 
sive coupling with coupling constant /c@: 

Rossler Lorenz (diffusively coupled) 

x\ = -(x 2 +x 3 ) y\ = a(y 2 - yx) 

x 2 = xi+ ax 2 y 2 = ryi - y 2 - 2/12/3 + k{x 2 ™ 2/2) 

x 3 = b + x 3 (xi - c) y 3 = yiy 2 -by 3 



V. CONCLUSION 



To conclude, we emphasis that the numerical study 
made in this paper encourages to take up further re- 
search on the apparent connection between BD and GS 
by divulging the facts how BDs can assist in assessing 
the robustness of GS and predictability of nature of 
the systems in GS. The classification of various chaotic 
systems in equivalence classes have been done. BD has 
been shown to be a distinguishing characteristic of such 
an equivalent class. In future, similar studies on chaotic 
maps and bi-directionally coupled synchronised chaotic 
systems would be of interest. 
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